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Abstract 

Higgs-boson production in association with bottom quarks, pp/pp bbH + X, is one of the most 
important discovery channels for supersymmetric Higgs particles at the Tevatron and the LHC. We 
have calculated the next-to-leading order QCD corrections to the parton processes qq, gg — > bbH 
and present results for total cross sections and for distributions in the transverse momenta of the 
bottom quarks. The QCD corrections reduce the renormalization and factorization scale depen- 
dence and thus stabilize the theoretical predictions, especially when the Higgs boson is produced in 
association with high-p^^ bottom quarks. The next-to-leading order predictions for the total cross 
section are in reasonable numerical agreement with calculations based on bottom-quark fusion 
bb^H. 



The Higgs mechanism P| is a cornerstone of the Standard Model (SM) and its supersymmet- 
ric extensions. The masses of the fundamental particles, electroweak gauge bosons, leptons, 
and quarks, are generated by interactions with Higgs fields. The search for Higgs bosons is 
thus one of the most important endeavours in high-energy physics and is being pursued at 
the upgraded proton-antiproton collider Tevatron with a centre-of-mass (_CM) energy of 
1.96 TeV, followed in the near future by the proton-proton collider LHC [3| with 14 TeV 
CM energy. 

Various channels can be exploited to search for Higgs bosons at hadron colliders. Higgs 



radiation off bottom quarks 

pp/pp ^ bbH + X via qq,gg-^bbH (1) 

is the dominant Higgs-boson production mechanism in supersymmetric theories at large 
tan/3, where the bottom-Higgs Yukawa coupling is strongly enhanced. The parameter 
tan/? = V2/V1 is the ratio of the vacuum expectation values of the two Higgs fields gen- 
erating the masses of up- and down-type particles in supersymmetric extensions of the SM. 
H = HsM, h^i may denote the SM Higgs boson or any of the CP-even neutral Higgs 
bosons of supersymmetric theories. In the SM the Yukawa coupling strength is given by 
gffH,SM = mf/v, where rrif denotes the mass of fermion / and v ~ 246 GeV is the vacuum 
expectation value of the SM Higgs field. In the minimal supersymmetric extension of the 
SM (MSSM) the Yukawa couplings are modified by the factors gfjn = 9 f fH, mssm/ 9 ffH, sm, 

cos a ^ sin a „ sin a ^ cos a , . 

9uuho = — — ^, guuHO = - — ^, gddho = ddciHO = ^, (2) 

sm p sm p cos p cos p 



where v = y f ^ + f I and a is the mixing angle of the CP-even Higgs fields and j5| . 

Leading-order (LO) predictions for the cross sections ((T)) are plagued by considerable 
uncertainties due to the strong dependence on the renormalization and factorization scales, 
introduced by the QCD coupling, the parton densities, and the running 6-quark mass in the 
bottom-Higgs Yukawa coupling. In this letter we present the first calculation of the next- 
to-leading order (NLO) QCD corrections to associated bbH production at hadron colliders 
through the parton processes qq, gg — > bbH. 

Note that the inclusive cross-section prediction for the process gg bbH develops poten- 
tially large logarithmic terms oc ln(mf,/Q). The logarithms arise from the splitting of gluons 
into bb pairs where the 6-quark mass acts as a cutoff for the coUinear singularity. The large 
scale Q corresponds to the upper limit of the coUinear region up to which factorization is 
valid. It has been argued that Q is of the order Mh/A or less P|. The hi{mh/Q) terms can 
be summed to all orders in perturbation theory by introducing bottom parton densities , 
provided that the 6-quarks are produced predominantly at small transverse momentum. In 
this approach [1, |^ [l3], the counting of perturbative orders involves both the strong cou- 
pling Q^s and the logarithm ln(mb/Q). The LO process is Higgs-boson production through 
bb fusion, bb H . The first order corrections comprise the 0{as) corrections to bb ^ H 
and the LO process gb bH, which is suppressed with respect to bb ^ H by 1/ ln{mb/Q). 
The second-order corrections, involving bb H at NNLO, have been completed recently 
[To] and strongly reduce the spurious scale dependence of the total cross-section prediction. 

Higgs searches in the bbH channel, however, do in general employ cuts on the 6-quark 
transverse momenta and thus require theoretical predictions for exclusive final states. This 
problem cannot be solved within the 6-quark density approach, since per definition the 



2 



transverse momenta of spectator 6-quarks are integrated out. The solution in this case is 
provided by calculations based on the parton processes gg, gg bbH. In this letter we 
present NLO results for both inclusive and exclusive observables, based on these parton 
processes. We also include a first numerical comparison of our results with calculations 
based on the bb fusion process. 

Generic diagrams that contribute to the processes ([[]) at LO are displayed in Fig. ^a). 
The NLO corrections comprise virtual one- loop diagrams, Fig. ^b), gluon radiation pro- 
cesses, Fig. life), and gluon- (anti) quark scattering reactions. Fig. [Tfd). The latter two add 
incoherently to the virtual corrections. The calculation of the NLO QCD corrections to 
the processes qq,QQ^ QQH, where Q denotes a generic heavy quark, has been described 
in detail in Ref. (see also Ref. [13 )• Iii the following we define the physical input pa- 
rameters and quote their values as chosen for the numerical analysis. The renormalization 
of the strong coupling as(/i) and the factorization of initial-state collinear singularities are 
performed in the MS scheme. For the calculation of the pp and pp cross sections we have 
adopted the CTEQ6L1 and CTEQ6M 3| parton distribution functions at LO and NLO, 
corresponding to = 165 MeV and Af^ = 226 MeV at the one- and two-loop level of the 
strong coupling as{fi), respectively.^ The top quark is decoupled from the running of as(/x). 
We evaluate the bottom-Higgs Yukawa coupling with the running 6-quark mass rnf,{n) de- 
fined in the MS scheme in order to sum large logarithmic corrections oc ln(mb/Mj:f) related 
to the renormalization of the Yukawa coupling [1^. The 6-quark pole mass has been set 
to mh = 4.62 GeV, corresponding to a MS mass frihijni,) = 4.28 GeV The strength of 
the SM Yukawa coupling is fixed by gbbH = 7nb{l^)/v, where v = {^/2Gf)~^^'^ is the vacuum 
expectation value of the Higgs field and Gp = 1.16639 x 10^^ GeV^^. 

In the following we present the NLO QCD predictions for associated bbH production at 
the Tevatron and the LHC in the SM. In LO, the transition to the MSSM can be done by 
a simple rescaling of cross sections and distributions with a factor g^^^^o /bbH" ' since the LO 
matrix element is proportional to the bottom-quark Yukawa coupling gbbho/bbHO- In NLO, 
this rescaling is spoiled by one-loop diagrams in which the Higgs boson couples to a closed 
top-quark loop (the corresponding quark loops of the first two generations are negligible). 
However, our calculation reveals that these contributions only amount to about — 5%(— 10%) 
of the NLO cross section at the Tevatron (LHC) in the SM. Thus, to a very good approxi- 
mation, also the cross sections and distributions including NLO QCD corrections follow the 
simple LO rescaling. We do not consider supersymmetric-QCD and electroweak corrections 
in this paper. ^ For a consistent comparison of LO and NLO results, all quantities {as, the 
parton densities, the running 6-quark mass, and the partonic cross sections) are calculated 
in LO and NLO, respectively. We first discuss the renormalization and factorization scale 
dependence. All scales (including the scale of the running 6-quark mass) have been set equal, 
H = fj,F = fiR, and are varied around the central value /iq = (2mb + Mh)/2. Figures El and 
01 show the scale dependence for both the total cross section and the cross section with two 



^ The CTEQ6 parton distribution functions involve five active quark flavours. In our calculation we do 
not include b quarks as active partons. The gluon distribution is enhanced in a scheme with four active 
flavours, but the effect is numerically small, less than 5%, and has been neglected in the present analysis. 

^ For large tan/3, supersymmetric loop corrections to the hhhP /H^ vertices may become important |l6| . In 
the present QCD analysis, the major part of these corrections can be absorbed into effective Higgs Yukawa 
couplings, since emission and reabsorption of virtual heavy supersymmetric particles is confined to small 
space-time regions compared with QCD subprocesses involving massless gluons. 
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TABLE I: Total cross sections for gg, gg — > bbH + X and hh ^ H + X [lo| at the Tevatron and 
the LHC. In the qq,gg bbH + X calculation the renormalization and factorization scales have 
been set equal to fj, = ij.o/2 = (2m;, + Mh)/^ and are varied between /i = /xo/4 (upper limit) and 
/i = /io (lower limit). In the calculation based on the bb ^ H + X process the renormalization scale 
has been fixed to iir = Mh, while the factorization scale is varied between O.IMh < fJ'F ^ 0-7 Mh 
with hf = 0.25Mff as the central scale. Note that the bb ^ H + X calculation employs MRST 
LO and NNLO parton distribution functions fl7|. 
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high-pr bottom quarks. The reduction of the spurious scale dependence at NLO is particu- 
larly striking for the exclusive cross section where both 6-quarks are required to be produced 
with pt > 20 GeV. A significant reduction of the scale dependence is also observed when 
demanding only one of the 6-quarks to be produced at large pt (not shown). The residual 
scale dependence of the inclusive cross section is still sizeable at NLO, and further work is 
needed to improve the theoretical prediction. Note that large logarithmic corrections may 
spoil the convergence of perturbation theory if /i is chosen too low. This is evident from 
the exclusive cross section at the Tevatron, Fig. |2l where the NLO prediction would even 
turn negative for /x < /xq/S. A similar observation has been made for tiH production at 
the Tevatron ^J]. We note that the bbH cross section at the Tevatron and the LHC is 
completely dominated by the gluon-gluon fusion and gluon-(anti) quark scattering reactions. 

The total cross section for associated bbH production at the Tevatron and the LHC is 
displayed in Fig. HI as a function of the Higgs-boson mass. We have set all scales equal to 
H = Ho/2 = {2mb + Mh)/4:, which is a suitable factorization scale choice for calculations 
based on the bb fusion process jESIlSl- ^'^^ ^^is choice of scales, the NLO QCD corrections 
increase the theoretical prediction significantly, by a factor 1.6 — 2.3 at the Tevatron and 
1.1 — 1.8 at the LHC. Representative results for the total cross section are listed in Table H] 
and are compared with the NNLO calculation based on 6-quark fusion bb ^ H y^]. When 
including higher-order QCD corrections, the results obtained in the two approaches are 
in reasonable numerical agreement.^ In particular, the cross-section predictions for bbH 
production at the LHC are compatible within their respective scale uncertainties. We would 
like to emphasize, however, that in spite of the qualitative numerical agreement between 
the qq, gg bbH and bb ^ H cross-section predictions, the intrinsic uncertainties of both 



^ Note that contributions involving Higgs radiation off a closed top-quark loop have not been included in 
the 6-quark fusion calculation 10]. While these contributions are negligible in supersymmetric theories 
at large tan /?, they do affect the SM cross section. Higgs radiation of top-quark loops has been included 
in the calculation presented in this letter and reduces the cross-section prediction in the SM by 5 — 10%. 
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approaches have not yet been fully quantified. 

In Figures |31 and IHl we show the cross section with one or two high-p^ bottom quarks 
as a function of the minimal 6-quark transverse momentum. We observe sizeable positive 
corrections if the 6-quarks are produced at small px- If both 6-quarks are required to be 
produced with larger > 20 GeV the NLO corrections reduce the LO prediction. The cross- 
section predictions with a single high-gr bottom quark will be compared with calculations 
based on the parton process gb bH j9| in a forthcoming publication. 

In summary, the strong scale dependence of the LO predictions for the processes pp/pp ^ 
bbH+X, which provide important search channels for supersymmetric Higgs bosons, requires 
the inclusion of QCD corrections. For the total cross section, we find a reduction of the 
spurious scale dependence at NLO, but further improvements are needed to quantify the 
accuracy of the prediction. Still, our results are in reasonable numerical agreement with 
calculations based on the bottom-quark fusion process. We find a drastic reduction of the 
scale dependence at NLO when the Higgs boson is produced in association with high-p^^ 
bottom quarks. The improved NLO predictions for the exclusive cross sections can thus be 
taken as a base for experimental analyses at the Tevatron and the LHC 



Note added in proof 

After submission of this paper, an independent calculation of the exclusive bbH cross 
section with two high-pj- 6-quarks has been presented in Ref. 3]. The authors of Ref. ^ 



have chosen a different set of input parameters and final-state cuts, so that a direct com- 
parison with the results presented in our paper is not possible. A systematic comparison of 
the two calculations is, however, now in progress, and good agreement has been found for 
the exclusive bbH cross section with high-p^ ^-quarks, see Ref. Ref. [l^ also includes 
a comparison of the cross-section predictions with calculations based on the parton process 
gb bH. 
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FIG. 1: A generic set of diagrams (a) for the Born level, (b) for virtual gluon exchange, (c) gluon 
radiation and (d) gluon- (anti) quark scattering in the subprocesses qq,gg bbH etc. 
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FIG. 2: Variation of the LO and NLO cross sections with the renormahzation and factorization 
scales for bbH + X at the Tevatron. 




FIG. 3: Variation of the LO and NLO cross sections with the renormahzation and factorization 
scales for pp —> bbH + X at the LHC. 
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FIG. 4: Total cross section for pp — > bbH + X at the Tevatron (upper figure) and pp — > bbH + X 
at the LHC (lower figure) as a function of the Higgs boson mass. 
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FIG. 5: Cross section for pp bbH + X at the Tevatron with one or two high-pr bottom quarks 
as a function of the minimal 6-quark transverse momentum. 
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FIG. 6: Cross section for pp — > bbH + X ai the LHC with one or two high-pr bottom quarks as a 
function of the minimal 6-quark transverse momentum. 
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